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We study the roton-like dip in the magnon dispersion at the boundary of the Brillouin zone in the 
isotropic S = 1/2 Heisenberg quantum antiferromagnet. This high-energy feature is sometimes seen 
as indication of a fractionalization of the magnons to spinous. In this article, we provide evidence 
that the description of the dip in terms of magnons can be improved significantly by applying more 
advanced evaluation schemes. In particular, we illustrate the usefulness of the application of the 
principle of minimal sensitivity in varied perturbation theory. Thereby, we provide an example for 
the application of this approach to an extended condensed matter problem governed by correlations 
which can trigger analogous investigations for many other systems. 



PACS numbers: 75.10.Jm, 75.30.Ds, 02.30.Mv, 75.50.Ee 



I. INTRODUCTION 

Quantum antifcrromagnctism is a long-standing issue 
which cannot be reviewed briefly, but see for instance 
Ref. 1. Yet there are still important aspects which are 
not clear. In particular, the dynamics at high energies 
is not yet fully understood. But the quantitative un- 
derstanding of the high-energy dynamics is of increas- 
ing experimental relevance. For instance, the dispersion 
of spin waves in the two-dimensional parent compounds 
of the high-temperature superconductors displays dips 
which can only be accounted for by considering subdom- 
inant exchange couplings as well, for references and fur- 
ther discussion see Ref. 2. The relevance of the magnetic 
dynamics in the total Brillouin zone including the high- 
energy behavior close to the Brillouin zone boundary has 
recently been emphasized by inelastic X-ray scattering 
for doped and undoped cuprates'^. 

The purpose of the present article is partly method- 
ological. Thus we discuss the fundamental system, 
namely the isotropic Heisenberg quantum antiferromag- 
net with spin 5 = 1/2 and nearest neighbor exchange 
J > 



H 



^1+6, 



(1) 



where the subscript i runs over all sites and S runs over 
the vectors to the adjacent sites. This model is very 
well studied. But we focus on the roton-like dip at the 
wave vector k = (tt, 0) and its equivalent values (the 
lattice spacing is set to unity). This feature represents 
an open issue because it eludes precise calculation within 
spin wave theory. 

On one hand, the dispersion w(k)/5 at k = (7r,0) and 
at k = (7i'/2, 7r/2) takes precisely the same value in linear 
spin wave theory which represents the leading order in 



an expansion in 1/S, where S is the spin value. In the 
next-leading, first order this remains true as well. On 
the other hand, series expansion for S" = 1/2 around the 
Ising limit predicts a dip of 8.6%'*'^ which is confirmed 
by quantum Monte Carlo with a dip of 9.6%^. 

The idea suggests itself that further corrrcctions of spin 
wave theory in 1/5* cure the above discrepancy^"^. But 
this seems not to be the case. Although in second order 
l/S'^, a small dip appears, it takes only 1.4% which is 
far from what one would like to have^''. This number 
improves by about a factor of 2 upon passing to the third 
order 1/5'^ to 3.2%. But this is still far from the series 
and quantum Monte Carlo result The convergence 
turns out to be particularly slow. 

One may view this observation as a mere mathematical 
problem. But one may also wonder why the convergence 
is so slow and come to the conclusion that the underlying 
physical description in terms of spin waves, also called 
magnons, is not appropriate and that the true nature of 
the elementary excitations is a different one, for instance 
that the magnons disintegrate to spinous. Such a view is 
indeed discussed in the interpretation of the experimental 
findings^^"^^. In Ref. 13 the dip is given to be 7(1)% 
analysing the experimental data. The naive analysis of 
the peak positions in Figs. 3a and 3b in Ref. 13 suggests 
even about 10% for the dip. In any case, a sizable dip is 
an experimentally well-supported fact. 

Concerning the issue of the elementary excitations, it 
is useful to recall a well-studied system where the same 
question was discussed. In Heisenberg S = 1/2 lad- 
ders with two legs the elementary excitations are S = 1 
triplons because no long-range order occurs. But multi- 
particle continua are sizeable as well^^'^^ and they may 
be taken as precursors of a fractionalization towards 
spinons^^. Interestingly, the important multi-particle 
continua and their energetic vicinity to the dispersion of 
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the elementary triplons""^^ induces a dip in their dispersion 
at fc = compared to fc = 7r/2. But high order pertur- 
bative results are necessary to capture this effect^^'^^'^^. 

These observations led us to look for a quantitative de- 
scription of the dip in the dispersion on the square lattice 
in terms of magnons. Since standard perturbation the- 
ory seems to be not particularly efficient, see above, we 
follow a modified approach. The basic idea is to stick es- 
sentially to a second order perturbative approach, but to 
vary the starting point of the perturbation. This means 
that we vary the unperturbed Hamiltonian Hq in order 
to improve the results. Of course, the total Hamiltonian 
H = Hq + Hp may not be changed so that a variation of 
Hq will automatically imply a variation of the perturbing 
part Hp. 

Let us assume that Hq depends on some parameter u 
or on a set of parameters u which we may vary. How docs 
one determine the appropriate starting point Ho{uo) for 
the perturbation? In this question we follow the prin- 
ciple of minimal sensitivity^^. The underlying idea is 
that the exact diagonalization of H docs not depend on 
the starting point uq. This, of course, will not be true 
for a generic approximative scheme. Thus a quantity 
such as the ground state energy Eq will depend on u if 
computed approximately: Eq^^'^Iu). Then, the principle 
of minimal sensitivity suggests to choose uq such that 
Eq^^'^^u) depends minimally on u in the vicinity of uq. 
Thus, one should choose local extrema or saddle points 
to determine the starting value tig. Given that E^^^\u) 
is differentiable one obtains as defining equation 



II. MODEL AND METHOD 
A. Model 

Expressed in the usual spin operators the Hamiltonian 
is given in Eq. (1). Exploiting that the square lattice is 
bipartite we may write 



H = J 



E 



(3) 



where Ta is the lattice made only from all A sites. Fi- 
nally, we will focus on 5* = 1/2. But for introducing the 
bosonic representation it is convenient to treat general 
spin. We use the Dyson-Maleev representation 



^Ai 



'2S 



a aiUi 



S 



Ai 



Shj = 

,(t) 



S 



a] a. 



'2S 

-5^ 



{2S) 



S 



Ai 



[23) 



S 



Bj 



^2Sa] 



'2Shi 



(4a) 



(4b) 



where a\ ' '' are bosonic creation/annihilation operators on 

the A-sites and ll'p on the S-sites. Next, we transform 
these bosonic operators in momentum space. We stress 
that the momenta k are taken from the magnetic Bril- 
louin zone (MBZ), which is a tilted square in fc-space with 
the corners (±7r,0) and (0,±7r), because the real space 
coordinate i runs over T a- 



0. 



(2) 



B. Basic Steps 



Note that for n parameters Ui the above prescription im- 
plies n equations dmE^^'^'' {u)\^^z^^ = to determine ug. 

If the perturbation is performed around a bilinear 
bosonic Hamiltonian the dependence Hq{u) may equiva- 
lently be replaced by the dependence of the set {ai(u)} 
of annihilation operators and the corresponding creation 
operators which diagonalize Hq{u). This approach has 
been used to illustrate the usefulness of the principle of 
minimal sensitivity in perturbative calculations and for 
continuous unitary transformations'^'^". Below, we will 
use it to improve perturbative spin wave calculations for 
the Heisenberg model on a square lattice. The spin op- 
erators will be represented as introduced by Dyson and 
Maleev so that the bosonic approach can be directly put 
to use. 

The article is set up as follows. In the following section 
II, we introduce the model and its bosonic representation. 
In particular, the variation of the bosonic description will 
be explained. In Sect. HI we present results for the vari- 
ation in two parameters. Results for the ground state 
energy and for the dispersion are shown. Finally, the 
article is concluded in Sect. IV. 



Next, we perform a conventional Bogoliubov transfor- 
mation respecting translational invariance 



a 



t _ 



= /k^k + ^-k = m-kttk + ^k/^-k, (5) 



where and 



3^ arc the new operators in which we 
express the Hamiltonian. The prefactors 1^ and mk can 
be chosen at will as long as they fulfil + = 1 where 
we assume them to be real. The freedom of choice for 
these prefactors provides us with the possibility to choose 
the starting point of the perturbation theory as described 
in the Introduction. 

We find it convenient to parametrize the prefactors by 



rl 



Mk 



a^k = 



2/-tk 
1 - Mk 



1 



Mk 



1/2 



1/2 



mk 



1 



Mk 



2mi< 



1/2 



-a-k^k 



(6) 



where /Xk can still be chosen freely as long as |/Xk| < 1 
holds. To elucidate the above parametrization we recall 
that the choice 



Mk 



1 



1 



7k 77(cos(fc2,) + cos(fcy)) (7) 
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leads to the standard linear spin wave description. We 
will come back to this point below. 

First, however, we express the Hamiltonian in the fields 



„(t) 



and A 



(t) 



(8) 



We split it according to 

H = Hci + H]ji + Hqi, 

where Hc\ = —AJS^N simply stands for the classical 
ground state energy; note that here is the number of 
^-sites. The second term iJbi stands for the part which 
stems from the bilinear terms if H is expressed in the 
original bosonic fields in (4). It reads 

-f^bl = £^01 + ^Dl + -f^Bl, 



where 



En 



8JS ^ llxk{xk - •yi^) 

keMBZ 



Hbi := 4:JS 



E 

keMBZ 



^ik(a[ak + /?k/3k) 



Hbi 4JS J2 5ik(4/3-k + h.c.). 

keMBZ 

The momentum dependent prcfactors arc given by 
^ik := - 2xi,ji, + xl) 
-Bik := liilk - Sxk + 7ka:;k)- 

Doing the same for the quartic part Hqi yields 



(9) 

(10a) 
(10b) 
(10c) 



(11a) 
(lib) 



with 



£^02 
A2 



-JNAl 



2 



E 

keMBZ 

E 

keMBZ 



^k(2;k7k 



fl"D2 := 2J V ^2k(4"k + /3^/3k) 



^^B2 



where we find 



2J ^2k(4/^L + h.c.), 

keMBZ 



A2k 
B2k 



A2 ■ Aik 
A2 ■ Bik. 



(13a) 
(13b) 

(13c) 

(13d) 



(14a) 
(14b) 



Of course, the simplicity of the last relation results from 
the simplicity of the original model which is characterized 
only by nearest neighbor couplings which are all renor- 
malized by the mean-field effects in the same way. 

The quadrilinear interaction part is given by the 
normal-ordered expression 



1234 



2n?i4«I«2/?-3«4 



+ 4Vl2La\a3(3Ul3^2 + 2V,%l3l^a3(3_2f3-i 

+ V/2l4/3-l/3-2a3«4 + l//2'34/3-4/3-3/3-2/3-l 



,(15) 



where the subscripts i = 1, 2, 3, 4 stand for the momenta 
hi and — i stands for — k^. The conservation of momen- 
tum in the lattice is ensured by the Kronecker symbol J^l 
which implies ki -|-k2 = ks -|-k4 modulo reciprocal lattice 
vectors from the reciprocal lattice of the A-sites, i.e., 
g G F^ means g = (jitt, m7r) with the integers n, m if 
the lattice constant of the original square lattice is set to 

unity. The vertex functions V\234 ''^^^ given explicitly in 
App. A. 

Now we can combine the diagonal parts in 



H 



D 



00 



Hul + Hd2 



(16a) 



00 



4 J 



^k(4«k + /3l/3k) (16b) 



keMBZ 

£^00 := H,i + Eoi+Eo2 (16c) 

= -4JS^N - J{4SA2 + Al)N (16d) 

wk = (5 + iA2)/^(l-2a-k7k + 4) (16e) 
and the perturbing part Hp in 

Hp := Hb+Hv (17a) 

Hb := 4J J2 5k(4/^lk + h.c.) (17b) 

keMBZ 

Sk := {S + ^A2)lUjk - 2xk + 7^x1) (17c) 



(12) where Hb in (17a) stems from the sum Hbi + ^^B2- 



C. Approximate Evaluation 

A straightforward procedure is to use standard pertur- 
bation theory in Hy, for instance in second order, to com- 
pute the ground state energy Eq and the dispersion w(k) 
in an approximate way. (Note the difference between Wk, 
the dispersion in the unperturbed Hamiltonian H]j and 
the dispersion CL'(k) of the full Hamiltonian.) First, we 
focus on the ground state energy because its stationar- 
ity (2) will determine {/Xk}. The correction AEb due 
to Hb can be easily computed to infinite order in Sk 
analytically by Bogoliubov transformation 



AEb = -2JN{2S + A2)(^2 - ^2), 



(18) 



where A2 is the value for A2 if we diagonalize the bilinear 
Hamiltonian from the very beginning, i.e., A2 = A2 as 
given by Eq. (13b) for /Xk = \/ 1 — 7k • The correction 
AEv involving Hy are much more complicated so that 
we determine them only in second order in Hy 



AEv 



S'i2illl2hh)'^ Vi2li ^43^21 



J 



N'^ S4 wi + a;2 + + UJ4 



(19) 



Formally, there could also be a second order contribu- 
tion which is linear in Hb and in Hy, but no such term 
contributes to the ground state energy. Hence the ap- 
proximate ground state energy Eq is given by 



Eq = Eqo + AEb + AE\, 



(20) 
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In the same fashion, we compute the dispersion. The 
influence of Hb is again taken into account in infinite 
order yielding 



uJB{k)=2J{2S + A2)^l--/l (21) 
The additional second order correction Sy(k) reads 
Sy(k) = SBy(k) + Ei/y(k) (22a) 



EBv(k) 
Syy(k) 



p 



N 

'Hi 



(3) 



(22b) 



2;;qs 
kp 



^ {IplqkfS'i 

p,q,s 

^^pqs ^qpk 

Wk + + + UJs 



kppk/ 

^^pqs ^qpk 

^k — Wp — Wq — ^^s 

(22c) 



so that the total approximate dispersion finally is given 
by 



;(k) =cjB(k) + Sy(k). 



(23) 



Note that in the quadratic correction ^^^(k) both per- 
turbing terms Hb and Hy enter. 



for u < |f 1/4 and 



/k 



(25b) 



ucos(fc:c) cos(fey) + u(cos(fc2;) + cos(fcj,) + 2) 



otherwise. The above choice is motivated by two ar- 
guments. First, we intend to include the cosine terms 
which go beyond nearest neighbor processes. The sim- 
plest choice are the two next-nearest neighbor processes 
included above. Second, /k may not become negative be- 
cause |/Ltk| can exceed unity in this case. This is partic- 
ularly important at the border of the magnetic Brillouin 
zone where 7k = 0. In addition, /k may not exceed unity 
because /ik should not change sign. This implies that u 
and V may not be chosen too large. 

On the boundary, i.e., for 7k = 0, we choose = q 
ky = TT — q and obtain for u < \v\/A 

fg = {v/2 - 2m)(1 - cos(2g)) + \v\ - v (26a) 

and u > |w|/4 

/q = (2u-t>/2)(l + cos(2g)). (26b) 

We see that /g > is ensured. Note for future reference 
that for u = v/A no dispersion along the boundary of the 
magnetic Brillouin zone occurs. 



D. Variation of Hd 

The equations similar to the above can be found in 
many previous approaches^'^"^". The main difference is 
that in the previous equations /Zk was chosen such that 
Hb vanishes or appeared only in subdominant orders in 
1/5. The equations above for arbitrary /Xk are more gen- 
eral. They allow us to vary what we call an ak or /3k ex- 
citation. Thereby, the diagonal part of the perturbation 
H]j is varied and we can apply the principle of minimal 
sensitivity by looking for stationary points of Eq as it 
results from the approximate calculation. 

Pursuing this line of argument we should vary /ik at 
each point in the magnetic Brillouin zone in the range 
1 ^ lMk|- This, however, is far too ambitious because 
of the macroscopic number of parameters to be varied. 
Thus, to simplify the approach we choose a particular 
parametrization of /ik which relies only on a small num- 
ber of parameters. In the present work, we want to illus- 
trate the approach in principle and restrict ourselves to 
two free para meters. Moreover, it is reasonable to choose 
/ik close to -\/l — 7^ which would correspond to the cor- 
rect solution in the limit 5* — ^ cx). Therefore, our choice 
is 

^k = (l-/k)\/l-7^, (24) 

where 

/k := (25a) 
V cos(fc^) cos(fcj,) + \v\ + u{cos{kx) + cos{ky) — 2) 



III. RESULTS 

The results presented below are evaluated for S = 1/2. 
First, we analyze the dependence of the ground state 
energy on the chosen parameters. It turns out that the 
most interesting parameter region is u, w > 0. 

A. Ground State Energy 

Fig. 1 shows the energy dependence of Eq as given by 
(20) on the two parameters u and v. Obviously, no dom- 
inant local minima or maxima catch our eye in the upper 
panel. In the lower panel, one can presume a saddle point 
in the center of the figure. 

Closer inspection of this range of u and v, see Fig. 2, 
shows that there is a line of small cusps given hy u = 
v/A as long as v is not too large, see also right panel of 
Fig. 3 below. In view of the definition of /k in (25) the 
appearance of such a cusp may not surprise. In addition, 
the line u = v/A is special since it makes any dispersion 
at the magnetic Brillouin zone boundary vanish. 

To elucidate the energy behavior more quantitatively 
Fig. 3 shows two perpendicular cuts through the energy 
landscape of Fig. 1. The left panel in Fig. 3 follows the 
line of cusps along u = v/A. Clearly, a local maximum 
appears which is located at wq = 0.2502(1). But the de- 
pendence through this point along a line perpendicular 
to u = w/4 clearly displays the cusp at which is a local 
minimum. Note, that the definition (25) is prone to yield 
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FIG. 2: (Color online) Cuts of Eo{u,v)/{2N) perpendicular 
to the line m = u/4 along the line (uo + Am, vQ — Ati/4). The 
cusps at Aw = for not too large vq are obvious. 
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FIG. 1: (Color online) Approximate ground state energy 
Eo{u,v) per site of the original lattice, i.e., Eo/(2N), as com- 
puted by Eq. (20). If the calculation were exact, Eo should 
be constant. Due to the approximations used this is not the 
case. The best strategy is to look for local extrema or saddle 
points because they represent points where Eq is stationary 
at least locally. Upper panel: Overall view, no extrema or 
saddle points are discernible. Lower panel: For (7 > 0, the 
energy landscape displays more structure and a saddle point 
can be presumed in the middle of the figure. Note that for 
clarity the color coding in the lower panel is different from 
the one in the upper panel. 



cusps as stated above. But it is a priori not clear that 
these cusps are extrema. Since the point at {vq/A,vq) is a 
local minimum in one direction, but a local maximum in 
the perpendicular direction we are not observing a local 
extremum, but a sort of saddle point. If Eq were differ- 
entiable, for instance if it were smeared out a tiny bit by 
convolution with a narrow Gaussian, it would display a 
usual saddle point very close to (wo/4, wq). We interprete 
the occurence of this special point on the line {v/A,v) 
as evidence that the optimum /Xk should not display a fi- 
nite dispersion on the boundary of the magnetic Brillouin 
zone, cf. Eq. (26). 

For the precise determination of vq, calculations are 
done for various system sizes with linear extensions L = 
24, 32, 36. The extrapolation of the position of the local 
maximum yields wq — 0.2502(1) and the energy value at 
this position is found to be Eq/{2N) = -0.66902(2) J. 
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FIG. 3: Cuts of Eo{u,v)/{2N) along special lines. Left panel: 
along the lines of cusps u = i;/4. Clearly a local maximum 
appears at vo = 0.2502. Right panel: This cut follows the 
line through the point (uo = i'o/4, uo), but perpendicular to 
u — w/4. It is given by (mq + Am, uO — Am/4). The cusp at 
Am = is obvious. 



These values should be compared to the quantum Monte 
Carlo result^i Eo/{2N) = -0.669437(5) J and to the sec- 
ond order result of a plain 1/S expansion* which reads 
Eo/{2N) = -0.66999J. (Note that this number is re- 
ferred to as "third order" in Rcf. 8 because the authors 
include the classical energy in their power counting.) If 
we take the Monte Carlo expansion as reliable reference 
the variation of second order perturbation theory could 
reduce the deviation from 0.0006 J to — 0.0004J which is 
a reduction by about 25%. To judge the improvement we 
point out that passing from rather simple first order per- 
turbative spin wave theory Eq / {2N) = —0.67042 J to sec- 
ond order Eq/{2N) = —0.66999 J improved the ground 
state energy by 44%. 

We also stress that the improved result Eq/{2N) = 
—0.66902(2) J is obtained by using equations of the same 
complexity as the second order equations. The add-on 
is just the variation of the unperturbed starting point. 
In higher orders, this variation becomes an even more 
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efficient tool, see Ref. 20 for tfie discussion of tlie example 
of the quartic oscillator. 

In the end, however, our goal is not to obtain esti- 
mates for the ground state energy in the first place. In 
the varied perturbation theory, the analysis of the ground 
state energy primarily serves the purpose to fix the un- 
perturbed starting point Hd- 



B. Magnon Dispersion 

Above, we have determined the optimum starting 
point Hd{u,v) by analysing the dependence of the ap- 
proximate ground state energy. We identified the op- 
timum starting point to be given by (vo/4, uq) with 
Wo = 0.2502(1) where an unusual saddle point appeared. 
Next, we use this starting point to analyze the magnon 
dispersion in general and the dip between the values at 
(tt, 0) and at (7r/2,7r/2) in particular. 
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FIG. 4: (Color online) Dispersion of the Heisenberg lat- 
tice model in (i) linear spin wave theory (LSWT) including 
first order corrections, (ii) second order perturbation theory 
around the LSWT solution, i.e., for (m, v) = (0, 0), (iii) in var- 
ied second order perturbation theory around {vo/4:,vo) with 
vo = 0.2502 which corresponds to the saddle point of the 
ground state energy. 

Fig. 4 depicts the corresponding result compared with 
results from first and second order perturbative spin wave 
theory. First, we find that the overall shape of all three 
curves is very similar. This was expected from the com- 
parison of various perturbative results, high order series 
expansion and quantum Monte Carlo data, see Fig. 1 in 
Ref. 10. 

Second, the dip at (tt, 0) relative to the dispersion at 
(7r/2,7r/2) is more pronounced in the varied perturba- 
tion theory. We find that the dip takes the relative value 
3.3(1)% which is rather precisely the value which Sy- 
romyatnikov found in the much more complex third or- 
der perturbation calculation^". It improves the second 



order result of 1.4% by more than a factor of 2 while it 
is still away from the about 9% of dip obtained by series 
expansion^ or quantum Monte Carlo^. 

In detail, we find a;(7r,0) = 2. 2881(1) J and 
w(7r/2,7r/2) = 2.3643(2) J. The latter value is very close 
to the series value a;'''=™"(7r/2, 7r/2) = 2. 385(1) J and 
to the quantum Monte Carlo value uj^^^'^{T:/2,n/2) = 
2.39 J. The former is still by about 5% too high compared 
to a;'''=™"(7r,0) = 2. 18(1) J and wQ'^^(7r,0) = 2.16J. So 
there is still some way to go. 

But to judge the numbers obtained by varied per- 
turbation theory we also compare to the ordinary sec- 
ond order perturbative numbers a;^"'^(7r,0) = 2. 3586 J 
and a;2"^(7r/2,7r/2) ^ 2.3920 J and to the third or- 
der perturbative numbers lu^'-''^ {tt , 0) — 2. 3241(2) J and 
a;3'''^(7r/2,7r/2) = 2.4007(2) J. Relative to these num- 
bers, the varied perturbative results represent an im- 
provement, in particular in comparison to the plain sec- 
ond order results. It must be kept in mind that the var- 
ied perturbation theory is based essentially on the same 
equations as the plain second order results. 

These findings show that the variation of perturba- 
tive calculations indeed reduces deviations to the exact 
results. In this way, improved results are accessible with- 
out resorting to more complex higher order calculations. 



IV. CONCLUSIONS 

In summary, we investigated the Heisenberg quantum 
antiferromagnet in terms of spin waves (magnons) based 
on the Dyson-Maleev representation. The primary goal 
was to determine the dispersion of the magnons. A sec- 
ondary goal was the determination of the ground state 
energy. 

The approach employed for evaluation is based on stan- 
dard perturbation theory. But we do not pursue a plain 
expansion in l/S. Instead, we choose the unperturbed 
Hamiltonian Hu, which serves as starting point, arbi- 
trarily among bilinear bosonic operators. In the present 
article, we did not exploit the full freedom of choice of 
such operators but investigated a parametrization with 
two variables (u, v) which remains close to the bilinear 
bosonic Hamiltonian of linear spin wave theory. Con- 
sidering more variables would have obscured the funda- 
mental principle of the approach and it would have been 
rather cumbersome on the technical level. 

Following the principle of minimal sensitivity, we 
search the parameter space (u, v) for stationary points, 
i.e., local extrema or saddle points, of Eq{u,v). Thus 
the approach is called varied perturbation theory. Such 
a point is indeed found at (uq = fo/4, wq) with vq = 
0.2502(1). Due to non-differentiability, it is not an ordi- 
nary saddle point, but a point with a cusp-like minimum 
in one direction and a differentiable maximum in the per- 
pendicular direction. 

At this saddle point, the value of the ground state en- 
ergy is closer to the exact value than the plain second 
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order perturbation. Furthermore, the magnon dispersion 
obtained at this saddle point displays a more significant 
dip of 3.3% which is again more than twice as large as 
the plain second order calculation provides. Only the 
plain third order calculation yields a comparable dip of 
3.2%. The agreement with high order series results and 
quantum Monte Carlo data is still unsatisfactory because 
these approaches find a dip of about 9%. 

We attribute the remaining discrepancy to the low or- 
der (here: second order) perturbative approach which we 
employed to calculate Eo{u, v). We expect that a varied 
third order approach enhances the dip to about 6 — 7% 
percent, combining the factors of 2 from the variation 
(this article) and from passing from second to third or- 
der (Ref. 10). 

Previous investigations of the simple model of a quartic 
oscillator have shown that the variation of the starting 
point combined with improved evaluation schemes such 
as higher order perturbation theory^^ or continuous uni- 
tary transformation^^ is capable to provide very good 
quantitative results. 

The progress achieved in this article is two-fold: On 
the methodological side, we introduced the principle of 
minimal sensitivity in the calculation for an examplary 
extended solid state system displaying important corre- 
lations. 

On the physical side, we provided evidence that the 
dip in the dispersion of the square lattice Heisenberg an- 
tiferromagnct with S" = 1/2 can be explained in terms 
of magnons if more advanced approaches are used. In 



our opinion, one does not have to resort to fractional- 
ization into spinons as sometimes discusscd^'^ in order to 
understand the dip. 

But we admit that a quantitative reproduction of the 
dip has not yet been achieved so that further work is 
called for. Promising improved approaches to reach this 
goal comprise third order perturbative approaches and 
continuous unitary transformations^^. 
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Appendix A: Vertex Functions 

The vertex functions are given by 



^1234 = XiX3-f{l ~ ^) + XiXij{l ~ 4) + X2X3-f{2 - 3) + X2X4:-f{2 - A) 

-xi"f{l) - X2"/{2) - xiX3X4^'y{l - 3 - 4) - X2X3X4"/{2 - 3 - 4), (Ala) 

(2) 

^1234 -2:37(2 - 3) - X4-f{2 - 4) - XiX2X3-f{l - 3) - XiX2X4-f{l - 4) 

+xiX2"f{l) + 7(2) + xiX2X3Xijil - 3 - 4) 4- 0:3x47(2 - 3 - 4), (Alb) 
^1234 = -2:17(1 - 3) - 2:27(2 - 3) - 2:12:3X47(1 - 4) - 2:22:32:47(2 - 4) 

2:12:37(1) -t- 2:22:37(2) + 2:12:47(1 - 3 - 4) + 2:22:47(2 - 3 - 4), (Ale) 
1x22:32:47(1 - 3) -I- .Ti.T27(l - 4) + X3.T47(2 - 3) + 7(2 - 4) 

X47(2) - xiX2X47(l) - .T37(2 - 3 - 4) - ^1X2X37(1 - 3 - 4), (Aid) 
X2X3X47(1 - 3) - xiX3X47(2 - 3) - xi7(2 - 4) - X27(l - 4) 
-xiX47(2) + .T2X47(1) + xiX37(2 - 3 - 4) + X2X37(1 - 3 - 4), 



V^^^ - xi 

'^1234 — 2-1 



^1234 



y^^^ - -X 

''1234 — 2 



y(7) 
''1234 



y(8) 
'^1234 



1x47(2) + .X2X47(1) + xiX37(2 - 3 - 4) + X2X37(1 - 3 - 4), (Ale) 
17(1 - 3) - X37(l - 4) - xiX2X37(2 - 4) - xiX2X47(2 - 3) 
-^7(1 - 3 - 4) + xiX27(2 - 3 - 4) + X3X47(1) + xiX2X3X47(2), (Alf) 

xiX47(l - 3) + xiX37(l - 4) + X2X37(2 - 4) + X2X47(2 - 3) 

-xiX3X47(l) - X2X3X47(2) - xi7(l - 3 - 4) - X27(2 - 3 - 4), (Alg) 
xiX47(2 - 4) + X2X47(1 - 4) + xiX37(2 - 3) -I- X2X37(1 - 3) 

-xi7(2) - X27(l) - xiX3X47(2 - 3 - 4) - X2X3X47(1 - 3 - 4), (Alh) 



VuL = 2:1x37(2 - 4) + X2X37(1 - 4) + xiX47(2 - 3) -f X2X47(1 - 3) 

-xi7(2 - 3 - 4) - X27(l - 3 - 4) - xiX3X47(2) - X2X3X47(1). 



(Ali) 
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